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ABSTRACT 

A  computer  program  was  developed  to  obtain  the  wind 
tunnel  wall  corrections  for  wing  angle  of  attack,  induced 
drag,  and  pitching  moment  in  incompressible  flow.   The 
vortex  lattice  method  is  used  for  computation  of  these 
correction  factors.   The  program  can  be  applied  to  wind 
tunnels  of  arbitrary  cross-sectional  shape,  and  wings  of 
any  desired  planform,  subject  to  the  constraint  of  straight 
leading  and  trailing  edges. 
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I.   INTRODUCTION 

In  the  process  of  wind-tunnel  testing  the  experimental 
procedure  itself  constrains  and  alters  the  flow  past  the 
wing  under  examination.   This  modified  flow  induces  errors 
which  must  be  corrected  if  accurate  results  are  to  be 
obtained.   The  major  interference  effects  in  subsonic  wind 
tunnels  result  from  the  constraining  influence  of  wind 
tunnel  walls,  disturbances  from  model  supports  and  measuring 
equipment,  and  turbulence  and  non-uniformity  in  the  airflow 
itself. 

Russell  [Ref.  10]  treats  the  problems  resulting  from 
model  support  interference  and  introduces  a  method  for 
reduction  of  this  effect.   The  differential  equations 
governing  flow  within  a  wind  tunnel  are  the  same  as  those 
in  free  air.   However,  the  presence  of  the  wind  tunnel  walls 
changes  the  outer  boundary  conditions.   The  result  is  a 
series  of  perturbations  of  the  velocity-potential  gradient 
normal  to  the  model  surface.   The  effect  is  pervasive  in 
that  almost  all  of  the  measured  aerodynamic  quantities  must 
be  corrected. 

Generally,  Prandtl  is  acknowledged  as  having  established 
the  foundations  for  research  on  wind-tunnel  wall  interfer- 
ence.  His  lifting-line  theory  and  concept  of  trailing  vor- 
tices opened  the  doors  to  much  experimental  investigation 
in  this  area.   The  resulting  studies  considered  two-dimensional 
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and  three-dimensional  effects,  open  and  closed  wind  tunnels, 
and  various  geometrically  shaped  test  sections. 

References  1  and  8  offer  a  comprehensive  treatment  of 
wind-tunnel  wall  interferences.   Much  tabular  and  graphical 
data  are  presented  for  different  tunnel  shapes  and  wing 
spans.   Pankhurst  and  Holder  included  much  data  on  octagonal 
wind  tunnels,  while  Pope  concentrated  more  on  a  uniformly 
loaded  wing. 

With  the  advent  of  sweptback,  slender  wings  designed 
for  high-speed  flight,  the  lifting  line  model  was  no  longer 
suitable.   In  these  situations,  the  wing  lifting  system 
may  be  represented  by  distributed  horseshoe  vortex  elements 
[Refs.  6  and  7].   In  Ref.  9,  data,  graphs,  and  formulas  are 
easembled  from  various  sources  and  systematically  analyzed 
for  various  types  of  wings  in  wind  tunnels.   Although  the 
mathematical  formulations  are  tedious,  the  resulting  graphs 
and  tables  summarized  very  well  the  data  obtained  from  each 
theoretical  treatment. 

Even  with  all  this  information  available,  it  is  still 
an  arduous  task  to  approach  tabulated  data  with  a  specific 
wind  tunnel  shape  and  geometric  configuration,  and  hope  for 
a  solution  to  the  wall  interference  problem.   It  is  the 
purpose  of  this  paper  to  develop  a  computerized,  closed- 
form  solution  to  this  problem.   The  design  goal  was  to 
develop  a  digital  computer  code  that  would  be  applicable  to 
wind  tunnels  of  arbitrary  cross-section  and  wings  with 
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varying  planforms.  Such  versatility  would  enhance  the 
accuracy  of  test  data  and  provide  a  valuable,  accurate 
tool  in  wind-tunnel  testing. 
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II.   METHODS  USED 

A.   GENERAL 

Wind  tunnels  provide  an  invaluable  tool  in  the  design 
and  performance  prediction  of  wings  and  wings-and-bodies . 
However,  the  wind-tunnel  testing  environment  introduces 
its  own  set  of  problems.   The  interference  action  from  the 
walls  of  the  wind  tunnel  alters  the  development  of  induced 
velocities  by  the  wing  lifting  system.   As  a  result,  downwash 
formation  is  suppressed  and  the  wing  aerodynamic  character- 
istics are  those  of  a  wing  having  greater  aspect  ratio.   The 
accurate  determination  of  wind  tunnel  wall  interference  is, 
therefore,  an  important  consideration  in  wing  experiments. 

The  classical  wall  correction  factor,  6,  for  the  case 
of  a  single  horseshoe  vortex  representing  the  complete  wing 
is  defined  by: 


Aa   =   5  -  CT  (1) 

c   L 


where 


A  -1  w 

Aa   =   tan   — 


Using  small  angle  approximations,  the  tangent  is  approxi- 
mately equal  to  the  angle.   The  wall-correction  factor  may 
be  expressed  in  terms  of  wing  circulation,  vortex  span  and 
average  induced  downwash  velocity  by: 
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a   -   w  c  /on 

6   "   2br"  (2) 


w 


Calculation  of  6  is  dependent  on  the  determination  of 
the  downwash  velocity  induced  by  the  wind  tunnel  walls. 
Classical  methods  involve  the  utilization  of  an  image 
system  of  vortices  which  cancel  the  effect  of  the  wing 
lifting  system  at  the  wall.   The  boundary  condition  of  no 
flow  through  the  wall  is,  therefore,  preserved.   This  system 
can  provide  an  accurate,  closed-form  solution  to  the  ground 
effect  problem.   However,  when  four  walls  are  present,  each 
wall  must  be  represented  by  an  image  system  which  affects 
other  image  systems,  as  well  as  the  model  itself.   The 
result  is  a  series  of  infinite  systems.   Truncation  of  this 
system  introduces  inaccuracies  of  unknown  magnitude  and 
significance. 

This  method  of  images  has  been  applied  to  tunnels  of 
rectangular  cross-section  with  success.   However,  it  is 
inadequate  when  arbitrary  tunnel  cross-sections  are 
encountered.   Add  to  this  a  wing  with  sweep  and  taper,  and 
the  problem  becomes  significantly  altered. 

Therefore,  a  new  method  was  sought  to  solve  this  problem 
specifically  for  the  NPS  3.5  X  5  ft.  octagonal  wind  tunnel. 
The  technique  chosen  was  presented  by  Joppa  in  Refs.  4  and 
5.   This  approach  lends  itself  very  well  to  computer  solution 
and  can  be  applied  to  wind  tunnels  of  any  cross-sectional 
shape. 
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In  this  study,  an  extension  of  Joppa ' s  vortex  lattice 
method  was  made  to  allow  the  treatment  of  wings  with  varying 
shape,  sweep,  and  aspect  ratio.   Such  treatment  requires  a 
modification  of  the  single  horseshoe  vortex  representation 
of  the  wing  lifting  system.   The  method  of  Weissinger  as 
presented  in  Ref.  6  provides  an  accurate  means  of  representing 
variable  planforms. 

B.   VORTEX  LATTICE  METHOD 

Since  the  method  of  images  was  considered  to  be  limited 
in  application,  it  was  necessary  to  find  another,  more  versa- 
tile method  of  representing  the  wind-tunnel  walls.   The 
problem  is  essentially  one  of  finding  a  vorticity  distribu- 
tion which  will  accurately  model  the  wind-tunnel  walls.   This 
vortex  system  must  satisfy  the  boundary  condition  of  no 
flow  through  the  walls. 

In  Ref.  4,  Joppa  presents  a  method  which  accurately 
models  the  wind-tunnel  walls.   This  was  accomplished  by 
representing  the  tunnel  walls  with  a  lattice  network  of 
rectangular  vortex  loops.   This  vortex  system  lies  in  the 
plane  of  the  tunnel  walls  and  satisfies  Helmholtz's  theorem 
that  a  vortex  filament  can  neither  begin  nor  end  in  the 
flow.   Fig.  1  shows  a  representation  of  this  system  as 
applied  to  the  NPS  3.5  by  5  ft.  wind  tunnel. 

Each  vortex  segment  in  a  given  rectangle  has  vortex 
strength  V . ,  and  has  a  corresponding  control  point  located 
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Figure  1.   Wind  Tunnel  Geometry 
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in  the  center  of  each  rectangle  where  the  no-flow  boundary 
condition  is  satisfied.   The  last  vortex  lattice  rectangles 
located  far  downstream  are  represented  by  horseshoe  vortices 
trailing  to  infinity.   The  reasoning  is  that  far  downstream, 
where  the  influence  upon  the  wall  vortex  lattice  array  of 
the  bound  vortex  elements  on  the  wing  is  small  relative  to 
the  trailing  vortex  elements,  one  finds  that  the  vortex 
strengths  of  the  lattice  terms  do  not  change  in  the  stream- 
wise  direction.   Consequently,  by  a  straightforward  cancella- 
tion argument,  one  can  visualize  replacing  the  vortex  lattice 
rectangles  by  horseshoe  vortices  trailing  to  infinity. 

The  mathematical  presentation  of  the  theory  is  based  on 
the  Biot-Savart  law  describing  the  velocity  induced  at  a 
point  by  a  vortex  segment.   Joppa  in  Refs.  4  and  5  follows 
the  mathematical  formulation  beginning  with  the  Biot-Savart 
law.   This  development  is  summarized  in  Appendix  B. 

The  general  matrix  representation  of  this  method  centers 
around  the  boundary  condition  of  no  flow  through  the  wall. 
As  mentioned  previously,  this  condition  is  satisfied  at  a 
control  point  located  at  the  center  of  each  vortex  loop. 
Therefore,  it  follows  that  the  velocity  induced  at  each 
control  point  by  the  vortex  loops  must  be  exactly  cancelled 
by  the  velocity  created  by  the  wing  lifting  system.   This 
is  presented  below  in  matrix  representation 


[ROL]  (r)   =   [WOL]  (T  }  (3) 

Li  W 
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which  allows  a  solution  of  vortex  lattice  strengths  in 
terms  of  the  wing  vortex  strengths. 


{r>   =   [ROL]"1  [WOL]  {T  }  (4) 

Li  w 


Note  that  the  [ROL]  matrix,  although  possibly  quite  large, 
is  square  and  hence  may  be  inverted. 

[ROL]  represents  the  wall  influence  coefficient  matrix. 
It  is  obtained  by  calculating  the  downwash  induced  at  each 
wall  control  point  by  a  unit  value  of  each  vortex  loop. 
Therefore,  ROL. .  is  the  downwash  at  wall  control  point  "i" 
induced  by  vortex  loop  "j".   [WOL]  is  the  influence  coeffi- 
cient matrix  describing  the  wing  effect  on  each  loop  con- 
trol point.   WOL. .  is  the  downwash  produced  at  control  point 
"i"  by  the  j    wing  horseshoe  vortex. 

By  inverting  [ROL]  and  premultiplying  with  [WOL] ,  the 
strength  of  the  wall  vortex  loops  can  be  found  in  terms 
of  the  wing  circulation.   The  inversion  of  [ROL]  can  be 
very  time-consuming,  since  it  may  be  a  very  large  matrix. 
However,  by  using  symmetry,  the  order  of  this  matrix  can 
be  reduced  by  a  factor  of  four.   This  requires  the  wing  to 
be  placed  on  the  horizontal  and  vertical  centerline  of  the 
wind  tunnel. 

C.   AIRFOIL  REPRESENTATION 

In  Ref.  5,  Joppa  stated  that  the  wing  can  be  repre- 
sented by  a  lifting  system  more  complicated  than  a  single 


20 


elementary  horseshoe  vortex.   In  order  to  model  swept  wings 
of  arbitrary  planform,  a  more  definitive  lifting  system 
model  was  employed.   Gray  and  Schenk,  Ref.  6,  present  a 
modified  Weissinger  method  designed  to  find  the  load  dis- 
tribution of  wings  with  arbitrary  planform.   The  wing  is 
represented  by  a  series  of  horseshoe  vortices.   The  midspan 
of  each  horseshoe  vortex  bound  element  is  located  at  the 
wing  quarter-chord  position.   The  condition  of  flow  tangency 
is  satisfied  at  the  wing  three-quarter-chord  point  in  the 
center  of  each  horseshoe  vortex.   Figure  2  illustrates  the 
configuration  chosen. 

Pearson  [Ref.  7]  modified  this  technique  by  allowing 
the  bound  vortices  to  be  swept,  and  including  chord-wise 
horseshoe  vortices.   For  this  presentation  it  was  felt  that 
unswept  bound  vortices  would  accurately  represent  the  lifting 
system  of  a  swept  wing.   In  addition  it  was  felt  that  little 
deterioration  of  accuracy  would  result  if  the  chordwise 
system  of  bound  vortices  was  combined  into  a  single  bound 
vortex  at  the  quarter-chord.   Figure  2  illustrates  the 
configuration  chosen. 

An  important  concept  in  this  representation  of  the  wing 
is  that  of  flow  tangency  at  the  three  quarter-chord.   The 
downwash  angle  at  the  wing  can,  therefore,  be  determined 
by  placing  a  control  point  at  this  position  in  the  center 
of  each  horseshoe  vortex.   The  rationale  behind  the  three- 
quarter  chord  flow  tangency  condition  is  a  product  of 
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Figure    2.      Wing   Geometry 
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two-dimensional  flow  considerations.   Extension  of  this 
technique  to  the  three-dimensional  wing  has  been  shown  to 
produce  experimentally  confirmed  results.   However,  it  is 
still  unclear  why  this  technique  works  so  well. 

The  wing  aerodynamic  influence  coefficient  matrix  was 
then  developed.   [AERO]  represents  the  downwash  induced  at 
each  0.75  chord  wing  control  point  by  a  unit  value  of  wing 
horseshoe  vortices.   The  matrix  formulation  for  induced 
downwash  formation  in  free  air  is  presented  below: 


(w>   =   [AERO]  {T  }  (5) 

w 


Dividing  the  above  expression  by  freestream  velocity,  one 
obtains  an  expression  for  0.75  chord  control  point  induced 
downwash  angles.   The  flow  tangency  requirement  equates 
this  column  matrix  to  the  wing  angle  of  attack  column 
vector. 


{£}       =   [AERO]  {b  (6) 

V  3c/4  V 


and 


(a>   =   {£> 


V  - 
3c/4 


i  2r 

It  should  be  noted  here  that  —  is  equal  to  — ,  from 
the  Kutta-Joukowski  Law.   Substitution  into  equation  6 
yields  the  following  results: 
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{£}     =   0.5  [AERO]  {£}  (7) 

3c  q  FA 

4 

Solving  equation  7  for  {— } ,  one  can  obtain  the  wing  span 
load  distribution.   Since  —  is  equivalent  to  the  sectional 

q 

value  of  C«c,  integration  of  the  —  distribution  vs.  span 
provides  a  direct  measure  of  the  wing  lift  coefficient,  C  , 

Li 

When  one  assumes  a  uniform  angle  of  attack  value  of  one 
radian  along  the  span,  the  resultant  integration  of  the 
span  load  solution  will  yield  the  wing  lift  curve  slope, 
CT  .   It  should  be  borne  in  mind  that  these  techniques  are 
predicated  upon  principles  from  thin  airfoil  theory  and 
in  the  limit  for  an  infinite  aspect  ratio  wing  will  yield 
a  wing  lift  curve  slope  of  2n-(rad   ). 


D.   DETERMINATION  OF  THE  WALL  CORRECTIONS 

The  classical  wall  interference  factor  6 ,  which  includes 
a  ratio  of  wing  to  tunnel  cross-section  area,  is  generally 
found  by  representing  the  wing  as  a  single  horseshoe  vortex. 
Since  the  present  formulation  utilizes  ten  horseshoe  vor- 
tices per  semi-span,  6  in  a  modification  can  best  be  found 
by  continuing  our  matrix  algebra. 


a-,   =   am  +  5C_ 
FA       T      L 


-FA   -   C^  (8) 

aFA 
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a   =  JL. 

T      CT 

T 


By  assuming  that  the  lift  coefficient  is  the  same  in 
both  tunnel  and  free  air,  the  following  equations  result: 


+  6  CT  (9) 


CL         CL 

aFA        aT 


Cancelling  the  lift  coefficient  and  solving  for  5: 


6   =   pr-i—  -   pi-  (10) 

FA         T 


The  free  air  lift  slope  was  found  in  the  previous  section. 
The  wing  lift  curve  slope  for  the  wing  in  the  wind  tunnel 
can  be  found  by  a  simple  extension  of  this  same  method. 
The  total  downwash  at  the  wing  control  points  is  equal  to 
that  induced  by  the  wing  added  to  that  created  by  the  tunnel 
walls.   The  effect  produced  by  the  tunnel  walls  is  repre- 
sented by  the  vortex  rectangle-on-wing  influence  coefficient 
matrix  [ROW] .   ROW. .  is  the  downwash  induced  at  the  i   wing 
control  point  by  the  j    wall  vortex  loop.   Multiplying  this 
matrix  by  the  vortex  loop  strength,  one  obtains  the  total 
wall  effect. 


{w>   =   [AERO]  {r  }  +  [ROW]  {r_}  (11) 

W  Lt 
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Substitution  of  equation  4  for  {r_}  forms  the  following 

Li 

set  of  equations: 


{w}   =   [AERO]  {T  }  +  [ROW]  [ROL]"1  [WOL]  {V    } 

w  w 


{-}     =   0.5  [[AERO]  +  [ROW] [ROL]  1 [WOL] ] {-}      (12) 
V  3c  *   T 


By  establishing  one  radian  angle  of  attack,  solving 
for  — ,  and  integrating  the  result,  the  wind-tunnel  lift 
curve  slope  can  be  determined.   It  should  be  noted  that  the 
correction  factor  obtained  from  this  process  is  different 
from  the  classical  5,  since  the  correction  factor  computed 
above  already  contains  the  wing  area  to  tunnel  area  ratio 
influence.   Therefore,  the  correction  to  angle  of  attack 
can  be  obtained  by  a  direct  application  of  this  6. 


Act   =   5  C_  (13) 

L 


In  most  literature,  the  correction  for  induced  drag 
is  obtained  by  application  of  the  same  5  as  that  used  for 
angle  of  attack. 


ACD  =   <5  CL2  (14) 

i 
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However,  in  this  presentation,  the  induced  drag  correc- 
tion is  calculated  in  a  different  manner  because  of  added 
available  aerodynamic  information.   The  spanwise  induced 
drag  distribution  in  free  air  is  given  by  the  following 
formula : 


{|>     =   [X^J    [DAERO]  {^}  (15) 


q  FA        q   FA  q 


FA 


[DAERO]  is  the  wing  aerodynamic  influence  coefficient 
matrix  computed  with  the  wing  control  point  now  located 
at  the  quarter-chord  position.   The  wall  effect  can  similarly 
be  found  by  including  additional  matrix  terms  for  the  wall 
effect  on  the  wing  quarter-chord  position. 


{1}    =   [v£j   [[DAERO]  +  [ROW]  [ROL]_1[WOL]  ]  {-}    (16) 

q  T     q  T  q  T 


Integration  of  equations  15  and  16  result  in  the  free 
air  and  wind  tunnel  induced  drag,  respectively.   The  induced 
drag  correction  factor  can  then  be  found  in  the  following 
manner. 


vA  =  K-c4 

FA 


\   -   Kt  < 
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K_.   =   C^    /C2 

1  (17) 

K     =   KFA  "  KT 


AC^   =   K  C2 
Di        LT 


When  testing  a  model  in  the  tail-on  configuration, 
the  presence  of  the  wind  tunnel  walls  serves  to  modify  the 
downwash  angle  induced  at  the  tail.   This  altered  flow 
field  can  have  a  significant  effect  on  the  pitching  moment 
coefficient.   The  following  increment  in  pitching  moment 
coefficient  is  generated  by  the  horizontal  tail: 


ACrnT   =   (^TL  Vc^  CL    >  aTL  (18) 

TL  w       aTL 


The  term  in  parentheses  is  approximately  equal  to  the 

3C 
stabilizer  effectiveness,  ^^  .   The  tail  angle  of  attack 

a  0  . 

1 

can  be  estimated  as  follows: 


aTL   =   a  "  £  +  5i  (19) 


At  a  given  value  of  wing  lift  coefficient  an  incremental 
change  in  downwash  angle  at  the  tail  due  to  the  wall  effect 
can  be  computed: 
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A£   =   £T  "  £FA  (20) 


Therefore,  the  correction  to  pitching  moment  coefficient  is 
given  by  the  following  relation: 


AC   =   C  .  Ae  (21) 

m      mi  v   ' 


w 
Since  the  downwash  angle  e  is  equal  to  ^   at  the  tail 

position,  e  for  free  air  and  for  the  wind  tunnel  can  be 

found  from  the  following  relations,  respectively: 


E__   =   0.5 [WOT]  {£} 
tA  q  FA 


em   =   0.5 [[WOT]  +  [ROT] [GAMAL] ]  {-}  (22) 

T  q  T 


[WOT]  is  the  wing  effect  on  the  tail  position  and 
[ROT] [GAMAL]  represents  the  wall  effect  on  the  tail. 
Equation  20  can  then  be  used  to  find  Ae.   The  pitching 
moment  correction,  as  evaluated  in  equation  21,  is  thus 
dependent  on  being  given  an  experimental  value  for  the 
tail  effectiveness  parameter. 
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III.   RESULTS  AND  DISCUSSION 

A.   COMPUTER  PROGRAM 

The  computer  program  developed  consists  of  two  major 
subroutines.   Subroutine  Loop  calculates  the  downwash  at 
any  point  due  to  a  vortex  rectangle.   Subroutine  HSHOE 
computes  the  downwash  at  any  arbitrary  point  due  to  a 
horseshoe  vortex.   HSHOE  is  used  for  both  the  wing  horseshoe 
and  those  of  the  wall  system  located  far  downstream  in  the 
tunnel  walls.   The  program  offers  two  basic  options.   One 
option  establishes  whether  ground  effect  or  wind-tunnel 
wall  corrections  are  to  be  calculated.   The  other  option 
indicates  whether  the  tail-on  or  tail-off  case  is  to  be 
considered. 

Each  computer  run  offers  the  inclusion  of  many  wings 
of  different  span,  aspect  ratio,  taper  ratio,  and  sweep. 
The  wing  parameters  are  read  as  input  data  after  the  inver- 
sion of  the  large  wall-on-wall  coefficient  matrix.   There- 
fore, with  a  given  tunnel  configuration,  a  variety  of 
different  wings  can  be  included  with  a  minimum  of  computer 
time  consumed. 

Integration  of  the  wing  span  loads  in  accomplished  by 
a  subroutine  obtained  from  the  NPS  source  library.   It 
should  be  noted  that  in  the  integration  scheme  the  assumption 
is  made  that  the  wing  loading  at  n  =  0.0  is  the  same  as  that 
for  n  =  .05b.   Any  inaccuracy  resulting  from  this  assumption 
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appears  to  be  negligible.   Appendix  C  describes  the  necessary 
input  to  the  program. 

B.   GROUND  EFFECT 

Originally,  ground  effect  calculations  were  conducted 
merely  to  check  the  accuracy  of  the  vortex  lattice  method. 
As  previously  stated,  the  method  of  images  provides  a 
closed-form  solution  to  the  ground  effect  problem.   There- 
fore, this  provided  an  excellent  opportunity  for  a  conver- 
gence check.   As  an  object  of  academic  interest  an  option 
for  ground  effect  computation  was  retained  in  the  final 
computer  program. 

The  results  of  the  convergence  check  showed  that  the 
vortex  lattice  method  could  very  accurately  compute  the 
velocity  induced  at  a  point  by  the  presence  of  a  wall.   By 
applying  the  theory  previously  developed,  the  lift  curve 
slopes  for  free  air  and  ground  effect  were  computed.   Results 
were  then  compared  with  those  presented  in  Etkin  [Ref.  3]. 
The  free  air  lift  curve  slope  was  found  to  correlate  extremely 
well  for  both  straight  and  sweptback  wings. 

Furthermore,  in  Ref.  3,  Etkin  defines  a  ground  effect 
ratio  as  the  ground  effect  lift  curve  slope  divided  by  that 
in  free  air.   The  results  are  tabulated  in  Table  1,  and 
Figure  3  compares  the  results  obtained  from  the  computer 
program  with  those  presented  in  Etkin.   As  can  be  seen, 
agreement  is  very  favorable. 
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Figure  3.   Ground  Effect  Ratio  Plotted  as  a  Function 
of  Wing  Height  Above  the  Ground 
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Several  computer  runs  were  conducted  at  different  wing 
heights  above  the  ground.  For  each  run,  sweep  was  varied 
from  0  to  40  degrees.  Results  are  tabulated  in  Figure  4. 
It  can  be  noted  that  the  ground  effect  ratio  initially 
increases  with  sweepback  but  returns  to  the  unswept  value 
at  about  30  degrees. 

As  a  note  of  caution,  care  should  be  exercised  when 
the  wing  gets  very  close  to  the  ground.   It  is  necessary 
to  greatly  decrease  the  vortex  loop  size  for  accurate 
results.   This  procedure  seems  to  provide  a  much  more  finely 
defined  ground  representation.   Otherwise,  calculations 
will  result  in  inordinately  high  values. 

C.   WIND  TUNNEL  APPLICATION 

A  convergence  check  of  the  wind  tunnel  wall  correction 
factor  was  obtained  by  applying  the  program  to  a  3  X  5  ft 

rectangular  tunnel.   The  correction  factor  obtained  was 

S 
converted  to  the  classical  5  by  including  the  factors  -^ 

and  57.3  (radians  to  degrees).   The  result  of  0.1216  was 
compared  to  that  presented  by  Joppa  in  Ref.  5.   The  wing 
geometric  configuration  presented  in  Ref.  5  consisted  of  a 
single  horseshoe  vortex.   The  value  of  0.127  from  Ref.  5 
correlates  very  well  with  the  results  of  this  study. 

Tables  2  and  3  present  the  data  obtained  from  compu- 
tations with  the  NPS  3.5  X  5.0  ft  octagonal  wind  tunnel. 
Computer  runs  were  conducted  using  wings  of  varying  span 
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and  sweep  angle.   Figures  4  and  5  illustrate  the  span  loads 
calculated  for  a  typical  wing.   it  is  interesting  to  note 
that  as  sweep  angle  increased,  negative  induced  drag  appears 
on  the  outboard  wing  stations.   Therefore,  total  induced 
drag  is  seen  to  steadily  decrease  with  increasing  sweep 
angle. 

As  sweep  angle  is  increased,  the  free  air  lift  curve 
slope  decreases.   Correspondingly,  the  corrections  to 
angle  of  attack  and  induced  drag  are  reduced.   Figure  6 
graphically  displays  the  data  for  a  straight  wing.   In  the 
past,  the  correction  factors  for  the  NPS  3.5  X  5.0  ft 
octagonal  wind  tunnel  were  obtained  by  assuming  an  elliptic 
cross-section.   The  correction  factors  presented  in  Ref.  1 
for  elliptic  cross-section  are  0.610  for  angle  of  attack 
and  0.010  6  for  induced  drag.   In  comparison,  the  vortex 
lattice  program  computes  values  of  0.7086  and  0.01075, 
respectively.   It  is  interesting  to  note  that  the  induced 
drag  correction  is  almost  identical,  whereas  the  correction 
to  angle  of  attack  is  slightly  higher  for  the  octagonal 
wind  tunnel.   Figure  7  illustrates  the  angle  of  attack 
correction  factor  with  varying  tunnel  span  ratios.   The 
correction  factor  can  be  seen  to  steadily  increase  with  the 
span  ratio. 

It  is  an  interesting  observation  that  the  value  of 
induced  drag  for  a  wing  span  of  5.0  ft  was  almost  zero. 
This  compares  very  well  with  airfoil  theory  in  that  a  wing 
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SPAN  RATION 


Figure  6.   Angle  of  Attack  Correction  Factor  Plotted 
as  a  Function  of  Span  Ratio 
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Figure  7.   Angle  of  Attack  Correction  Factor  Plotted 

As  A  Function  Of  Sweep  Angle 
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which  spans  the  wind  tunnel  represents  an  approximation  to 
the  two-dimensional  wing.   Induced  drag  for  a  two-dimensional 
wing  should  equal  zero.   In  addition,  the  wind-tunnel  lift 
curve  slope  is  within  0.0006  of  the  theoretical  two-dimensional 
value  of  2tt  rad 

A  sample  computer  run  was  made  for  the  moment  correction 
factor.   The  tail  quarter-chord  was  placed  two  feet  behind 
the  wing  bound  vortex,  on  the  tunnel  centerline.   The 
correction  factor  obtained  for  a  straight  wing  of  aspect 
ratio  6.0  was  -0.2297.   For  a  wing  with  sweep  angle  of 
30  degrees  and  taper  of  0.3,  a  value  of  -0.1887  was  calcu- 
lated.  The  above  correction  factors  must  be  multiplied 
by  the  stabilizer  effectiveness  parameter,  using  consistent 
angular  measures,  in  order  to  find  the  pitching  moment 
correction. 
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IV.   CONCLUSIONS  AND  RECOMMENDATIONS 

It  is  felt  that  the  vortex  lattice  computer  program 
as  presented  in  this  paper  offers  a  very  reliable  and 
flexible  tool  for  the  determination  of  wind-tunnel  wall 
corrections.   It  is  capable  of  treating  wings  of  any  span, 
taper  ratio,  aspect  ratio,  and  quarter-chord  sweep,  subject 
to  the  constraint  of  straight  leading  and  trailing  edges. 
The  geometric  shape  of  the  wind  tunnel  cross-section  can 
similarly  be  varied  at  the  discretion  of  the  program  user. 
Resultant  accuracy  appears  to  be  excellent.   However,  cau- 
tion should  be  exercised  in  determining  the  vortex  loop  size 
As  the  wing  is  positioned  closer  to  the  walls,  a  smaller, 
more  definitive  vortex  loop  representation  must  be  used. 

A  logical  extension  of  this  program  would  be  the  inclu- 
sion of  wings  with  an  antisymmetric  load  distribution. 
Furthermore,  slight  modifications  could  incorporate  the 
ability  to  consider  the  wall  effects  on  other  stability 
derivatives.   In  its  present  form,  this  program  offers 
the  wind  tunnel  user  a  versatile  tool  for  determination  of 
wind  tunnel  wall  corrections. 
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APPENDIX  A 

The  general  assumptions  used  in  the  compilation  of 
this  program  are  listed  below. 

1.  The  two-dimensional  lift-curve  slope  is  2tt  rad 

2.  The  wing  load  distribution  is  symmetrical. 

3.  Aeroelastic  effects  were  not  considered. 

4.  The  wing  is  located  on  the  horizontal  and  vertical 
centerline  of  the  wind  tunnel. 

5.  The  leading  and  trailing  edges  of  the  wing  are 
straight. 
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APPENDIX  B 
DEVELOPMENT  OF  INDUCED  VELOCITY  EQUATIONS 

The  equations  required  to  solve  for  induced  velocities 
consist  of  two  main  types.   One  equation  solves  for  the 
velocity  induced  by  an  arbitrarily  oriented  vortex  segment, 
whereas  the  other  basic  equation  applies  to  horseshoe  vortex 
elements .   These  equations  have  been  incorporated  in  sub- 
routines LOOP  and  HSHOE,  respectively.   The  development  of 
these  equations  presented  below  follows  closely  that  found 
in  Ref.  4  and  5. 

According  to  the  Biot-Savart  law,  the  velocity  induced 
at  a  point  p  is  found  from: 


p  

W   =   t-TT  (COS  B.   +  COS  Bn)   V  (IB) 

47rh       1        2 


For  vortex  segments  oriented  in  the  Y-Z  plane, 


R   +  R 
cos  Bl  +  cos  B2   =  jj-g-gi  (S2-  CR-l-V2)   (2B) 


It  is  now  necessary  to  consider  v,  the  unit  vector 
which  establishes  direction.   From  Fig.  8,  it  can  be  seen 

that 


v  = 


R,  X  S 

J: -—  (3B) 

|r1  x  s| 
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Figure  8.   Velocity  Induced  By  A  Vortex  Segment 
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Noting  from  Fig.  8  that 


hs   =   |  Rx  X  S  |  (4B) 


then 


hs   =   [(R,  S   -  R,  S  )2  +  (R.  S   -  R.,  S  )2 
ly  z     lz  y       lz  x    lx  z 

+  (R,  S   -  Rn  S  )2]1/2 
lx  y     ly  x 


Calculating  the  induced  velocity  in  the  direction  of 
outward  normal  of  point  P  per  unit  r, 


Y     =  j   •  n  (5B) 


where  the  unit  normal  vector  is  defined  by: 


n   =   oi  +  (sin  9p)j  +  (cos  9p)k 


Therefore,  one  finally  obtains: 


=      Rl    +    R2       [s2    _    (R        R   )2] 

87TR1R2(hs)  ± 


(6B) 


[(sin   9p)  (RlzSx-  RlxSz) '+  (cosep)  (RlxSy-  RlySx)  1 
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A  conditional  statement  must  be  included  in  the  program 
to  set  the  downwash  equal  to  zero  if  hs  is  equal  to  zero. 
This  will  preclude  having  a  zero  term  in  the  denominator. 

The  trailing  elements  of  a  vortex  loop  can  be  treated 
with  more  simplicity.   For  these  segments, 

<r2  -  h2)1/2 

cos  B. 


1  R- 


(R,2  -  h2)1/2 
oos  B2   =   51 


However,  it  can  be  seen  from  Fig.  8  that 


?  2  2  2 

R    =   (xp-xA)   +  (yA-yp)   +  (zA-zp) 


h2   =   (yA-  yp)2  +  (zA-  zp)  2 


Combining  equations  7  and  8, 

xp  -  xA 

cos  B,   =   — ^-5 

1       Rx 


xB  -  xp 
COS  B2   =   -R^ 


(7B) 


(8B) 


(9B) 


As  a  result,  the  equation  for  the  trailing  elements 
of  a  vortex  loop  is  simply 
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cos  B1  +  cos  B 

Y     =   2 [(sin  9p)  (zA  -  zp) 

4irh 

(10B) 

+  (cos  6p) (yp  -  yA) ] 


For  the  horseshoe  vortex,  the  equation  for  the  bound 
segment  is  identical  to  equation  6.   Similarly,  the 
trailing  elements  can  be  treated  as  in  equation  11  with  a 
slight  modification.   In  this  case,  the  cosine  of  one 
enclosed  angle  is  equal  to  1,  since  the  trailing  element 
extends  downstream  to  infinity.   The  induced  velocity 
equation  becomes 


1  +  cos  B, 

Y     =     , 2 -[(zA  -  zp)j  +  (yp  -  yA)k] 

4?rh 


w      1  +  cos  B, 

r  =  * -[(sin  6p)  (zA  -  zp) 

4iTh 

(11B) 

+  (cos  8p) (yp  -  yA) ] 
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APPENDIX  C 
COMPUTER  PROGRAM  INPUT 

Card  Number   Column  Number    Description  Format 

1  1-72        Program  title  18A4 


1- 

-5 

MM=Number  of  Spanwise 
Stations 

15 

6- 

-10 

NN=Number  of  Downstream 
Stations 

15 

11- 

-20 

Deltax=Loop  length 

F10.0 

1 

NOPT=l  for  Ground 
Effect 

11 

2 

NTAIL=1  for  tail-on 

11 

1- 

-72 

Tunnel  y-coordinates 
(one  quadrant  only) 

8F10.0 

1- 

-72 

Tunnel  z-coordinates 
(one  quadrant  only) 

8F10.0 

1- 

-10 

Span 

F10.0 

11- 

-20 

Aspect  Ratio 

F10.0 

21- 

-30 

Taper  Ratio 

F10.0 

31- 

-40 

c/4  Sweep  (degrees) 

F10.0 

1- 

-10 

x-Coordinate  of  tail 
quarter-chord 

F10.0 

11- 

-20 

z-Coordinate  of  tail 

F10.0 

quarter-chord 

Input  cards  5  and  6  can  be  repeated  for  different 
planforms  and  tail  locations. 
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